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Abstract 

We put forward the concatenation of Quantum Arnold Transformations as a tool to obtain the wave 
function of a particle subjected to a harmonic potential which is switched on and off successively. This 
simulates the capture and release process of an ion in a trap and provides a mathematical picture of 
this physical process. 
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1 Introduction 

In a recent paper [Tj the authors proposed a transformation that maps states and operators from 
quantum systems in one dimension described by generic, time-dependent, quadratic Hamiltonians to the 
system of the free particle. This transformation is the quantum version of the Arnold transform (QAT), 
which in its original classical version [2] maps solutions of a certain type of classical equation of motion, a 
non- homogeneous linear second order ordinary differential equation (LSODE), to solutions of the classical 
equation for the free particle. Through this transformation, the class of quantum systems with quadratic, 
time-dependent Hamiltonian can be connected and its solutions and relevant operators related by using 
the free particle system as an intermediate point. This way, the symmetry properties of the free particle 
system are transferred to this whole class of quantum systems. 

In fact, the QAT proves to be useful in the construction of wave functions, the quantum propagator 
or the evolution operator, taking advantage of the simple properties of the free particle quantum system. 
Even more, some insight in the physics of the free particle can also be obtained by importing to this system 
features that are, in principle, naturally attributed to the quantum harmonic oscillator. This is the case 
of coherent and squeezed states in one spatial dimension, as well as Hermite-Gauss and Laguerre-Gauss 
states in higher dimensions, as it was shown in [3]. 

There are some related strategies in the literature employed to solve the Schrodinger equation with a 
general quadratic, time-dependent Hamiltonian. We could mention the method of looking for integrals 
of motion showed, for instance, in [4] and [21[6]. Some other references may be found in [TJ[3]. 
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A particular case of time-dependent potential is that of a harmonic potential which is turned on 
and off successively, preserving or not the original frequency oj of the oscillator. This could be useful to 
simulate the processes of capturing or releasing ions in traps. Usually, ion traps, like the Penning of the 
Paul trap, have time varying, periodic frequencies (for instance, for the Paul trap the classical equation 
of motion is a Matieu equation), but we shall restrict ourselves to constant frequency to focus on the idea 
to turning on and off the potential without introducing unnecessary technicalities. The purpose of this 
paper is to take advantage of the QAT that relates the harmonic oscillator and the free particle to obtain 
a method to compute analytically the wave function and the evolution operator for such a system. 

The complete time evolution operator along a lapse of time in which a particle is released and captured 
several times in a harmonic potential will be a product of evolution operators for the free particle and 
the harmonic oscillator. This is a consequence of the sudden or diabatic approximation if the processes 
of switching on (and off) are fast enough. The harmonic-oscillator evolution operation on certain states 
(namely eigenstates of the harmonic oscillator Hamiltonian) may be simpler than the free one. It will 
be shown that it is possible to split each free time evolution operator in two operations: a harmonic 
evolution and a quantum Arnold transformation, in such a way that the complete evolution can be 
written in terms of harmonic evolutions and QATs only. This might not be always the easiest way to 
compute the resulting wave function. We believe nevertheless that it is interesting to have the possibility 
that we present here at our disposal. This technique has also the advantage of being easily generalizable 
to generic quadratic potentials with abrupt jumps in its time-dependence. However, we will keep the 
simplicity of the harmonic potential to keep the idea as neat as possible. 

The paper is organized as follows. In Sec. [2] we give an overview of the Quantum Arnold Transfor- 
mation for general quadratic Hamiltonians. In Sec. [3] we provide the particularization for the case of the 
harmonic oscillator and fix some notation. In Sec. |4] we obtain the alternative expression for the time 
evolution operator in terms of QATs. A summary and some comments are given in Sec. [5j 



2 The Quantum Arnold Transformation 



The quantum Arnold transformation (or the inverse) relates the Hilbert space Tit of solutions of the 
free Schrodinger equation 

i^— = - — — (1) 

dt 2m dx 2 

associated with the corresponding classical equation 

x = 0, (2) 
to that space 7i' t , of solutions of the Schrodinger equation 

where the quantum theory of a generic LSODE0 

x' + fx' + wV = (4) 



1 We consider here the case of homogeneous classical equation and the corresponding quantum theory. For the presence 
of an external force, see pQ. 



2 



is realized. The quantities / and uj are, in general, dependent on time t'. The classical equation can 
be derived from a Hamiltonian function 

H = ?-e-f + (-mio 2 x' 2 )ef , (5) 

and, according to the standard canonical prescriptions, it leads to the Schrodinger equation ([3]). For 
/ linear in time and constant uj this equation is commonly known as Caldirola-Kanai equation for the 
damped harmonic oscillator 

Being both spaces of solutions of ([T]) and ([3]), Ht and %L respectively, related, we will also have 
the basic, constant of motion quantum operators associated with the classical functions position and 
momentum connected and realized as well-defined, constant of motion operators on both H t and H' t ,. 
However, it is important to note that the Hamiltonian of one of the systems is not transformed into the 
Hamiltonian of the other: the generators of the symmetry of the free particle, the Schrodinger group, are 
mapped into a set of operators closing the same Lie algebra in the system corresponding to the generic 
LSODE. The Hamiltonian operator of the LSODE-system is not the same abstract element in this algebra 
or, even, it does not necessarily belongs to this symmetry algebra. 

This is a generalization of the classical Arnold transformation A, and is obtained by completing A 
with a change of the wave function. While the classical Arnold transformation A is given by: 

A: Rxf — >RxT 

(a-'.i')— ► ( x ,t)=A[(x',t')) = 

where T" and T are open intervals of the real line containing t' = and t = 0, respectively, and iti(i') and 
U2(t') are independent solutions of the classical equation of motion in (x',t') ([4]), the QAT is written: 

A : T-L^i — 7- l~Lt 

0'(a/;* )^ ¥>(z, t) = A(<f> (x , t')) = A* ( vM*') e 5 h W)»2V') X <P'(x',t')) , 

where A* denotes the pullback operation corresponding to A, and W(t') = U\U2 — u\U2 = . It is 
straightforward to check that by this transformation the Schrodinger equation of the free particle is 
transformed into ([3]) up to a multiplicative factor which depends on the particular choice of the classical 
solutions ui and U2 (partial derivatives must be changed by the classical part of the transformation while 
wave functions are shifted by the quantum part). 

We impose on u\ and U2 the condition that they preserve the identity of t and x, i.e., that (x,t) 
coincide with (x',t') at an initial point t' , arbitrarily taken to be t' = 0: 

ui(0) = 0, u 2 (0) = 1, «i(0) = 1, n 2 (0) = 0. (8) 

This fixes a unique form of A for a given "target" LSODE-type physical system. However, the quantum 
Arnold transformation is still valid if solutions u\ and U2 do not satisfy ([8]) (see [I] for details). 

It will be useful to have a pictorial representation of the situation. There will be a common Hilbert 
space H of wave functions, which plays the role of initial values for both the solutions (f>'(x', t') E %' t , and 
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(6) 



(7) 



<p(x,t) £ Hti evolution operators, U(t,0) = U(t) and U'(t',0) = U'(t'), acting on the initial wave func- 
tions; and the quantum Arnold transformation A that relates Schrodinger equations and basic operators: 



T~Lt 4 Hfi 



u(t) 



U>(t>) (9) 



Ho = T~L y H = T~Lq 

i 

where 1 is the identity operator. Had the solutions u\ and U2 not been chosen satisfying ([8]), the price 
would have been that the relation in the lower part of the diagram ([9]) above would no longer be the 
identity and basic position and momentum operators would then be mixed (see pQ). 

Let us stress that QAT can be useful to quickly perform some calculations, avoiding tedious, direct 
evaluations which can become extremely involved in the system under study. For example, it can be used 
to compute the quantum propagator for any LSODE-type quantum system, following the idea of Takagi 
in [9] for the simple case of the harmonic oscillator, or even the evolution operator U'(t'), which becomes 
very difficult to evaluate exactly when the Hamiltonian is time-dependent and does not commute with 
itself at different times. 

Actually, the evolution operator of a LSODE system can be related with the free evolution operator. 
Having in mind the diagram ([9|), we write: 

A{U'{t')4>\x')) =U(tMx). (10) 

Here (f>' and ip are the same function of only one argument (x or x') and we will denote ip = (/>' = ip. 
Then, 

1 % 777, l Ti2 /2 

U'{t!)i>(x') = A- l {U{t)^(x)) = -j=e^ w ^ x ^*- 1 (C/(t)V'(x)) = 

- — e -2-w- 2 * A*- 1 (;7(t))A*- 1 (^(x)) . (ii) 



'U2 

To factorize the function ip and single out the general action of U'(t'), we compute 

A*- 1 ty(aO) = V(£) = e lo ^ u ^'^(x') , (12) 

where ^°^ 1 / U2 "> X a^ 7 " is a dilation operator which is not unitary. To unitarize this operator, the generator 
must be shifted (symmetrized) from X ' to x'-^j + ^, so that the true unitary operator is then 

U D (±) = e iog(i/« 2 )(*'^+i) = _L e iog(iMK^ _ (13) 

U2 \/U2 



But the factor ^= is already present in the previous expression of U'{t'). Therefore, it now reads 
U'(t') = e ^^ x ' 2 A*~\lJ(t))tj D {i- 2 ) = 

1 e ff^^' 2 e ^^ 2 ^ e iog(iMK^. (14) 
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It is remarkable that we have been able to write an exact expression for the evolution operator as a 
product of operators. No perturbative approximation method, which could become cumbersome in some 
cases, is needed for any LSODE-related quantum system to obtain the evolution operator. 

3 A particular case: the harmonic oscillator and the free particle 



Now let us give the specific expressions for the harmonic oscillator, that is to say, / = and w(i') = w, 
constant. The two independent classical solutions can be chosen as 

ui(t') = — sin(ujt') 
^2(0 = cos(wt') 

with W(t') = 1. It can be checked that the change of variables in the classical Arnold transformation 
results in: 

t' = - tan' 1 (ut) 

id 

x = cos(tan _1 (wi))x 



VI + 0J 2 t 2 

The solutions u\ and ui satisfy the conditions ([8]). Therefore, the diagram ([9]) can be simplified to: 

A 

H « Who 




(15) 



where U and ?7ho stand for the evolution operators of the free particle and the harmonic oscillator, 
respectively, while %q is the Hilbert space, either for the free particle or the harmonic oscillator, of 
solutions of their respective Schrodinger equations at t = 0. 

It is possible to apply the QAT, for instance, to the time-dependent eigenstates of the harmonic 
oscillator Hamiltonian H^ , with energy E n = Huj(n + |) and belonging to %ho> 



(16) 
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where Af n = (^) 4 yp^j • We obtain the following set of states, solutions of the Schrodinger equation 
for the free particle: 

x 2 5* i 

* M -^ e ^Q"'^7m>- (17) 
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where, in order to obtain a more compact notation, we have introduced the quantities L = y , with 

dimensions of length, and r = = with dimensions of time. We also denote by 5 the complex, 
time dependent, dimensionless expression 5 = 1 + not = 1 + i = 1 + £7- We have also used the 
fact that e~ lu)t ' = e~ ltan ( ut ) = ™. Note that with these definitions the normalization factor M n can be 

written as A/;, = ^Ll. 

This set of states constitutes a basis for the space of solutions of the free Schrodinger equation, since 
it is mapped from a basis for the harmonic oscillator through A, which is unitary. The first state of this 
basis, the one mapped from the harmonic oscillator vacuum state, is given by: 

= ^ e"4lW = 1^=1 e -ik , (is) 

which is nothing other than a Gaussian wave packet with center at the origin and width L. The parameter 
r is the dispersion time of the Gaussian wave packet (see, for instance, [313]). The re-th state represents 
a wave packet with n + 1 humps (see [3] ) . 

The family of wave functions (|17j) has been known in the literature as Hermite-Gauss wave packets 
and they have been widely used, in their two dimensional version, in paraxial wave optics [12]. Note that, 
making use of the classical solutions only, and through the QAT, we have been able to import the time 
evolution from the stationary states of the harmonic oscillator, ip' n (x' , if), to the non-stationary ones, 
vp n (x,t), without solving the time-dependent Schrodinger equation. This will be the motivation for the 
next section. 



4 Switching on and off the harmonic potential 

Time evolution of an eigenstate of a given Hamiltonian, and in particular the harmonic oscillator, 
is encoded in just a time-dependent phase factor in the wave function. We have seen that if a physical 
system S, defined by another Hamiltonian, can be linked to the harmonic oscillator through an Arnold 
transformation, then the quantum time evolution of S can be found performing a change of variables 
in the wave function of the harmonic oscillator and multiplying by a specific phase and factor, that is, 
carrying out the QAT, provided the necessary classical solutions are known. This way, it is possible to 
avoid solving the Schrodinger equation or computing the evolution operator whenever the initial state 
of the system S is an eigenstate of the harmonic oscillator Hamiltonian. The same observation can be 
made if we have knowledge of the time evolution of a certain state of the harmonic oscillator and we use 
it as a initial state of S. To illustrate the following reasoning, we are going to ascribe to the free particle 
system the role of S and use the results of the previous Section. 

The physical situation we intend to describe is that of a particle in a harmonic potential which is 
switched off at a given time To. After that, the particle will evolve freely and then, at time T%, is captured 
again by another (in particular, the same) harmonic potential. Time evolution will be harmonic again 
until it is released once more at T2. Finally, it will be detected at time T3. 

This kind of process is of practical relevance, for instance in the preparation of the 1-dimensional 
Hermite-Gauss free states (|17p . This idea was proposed in [13], and was named "Quamtum Sling". The 
vacuum state of the harmonic oscillator, when switched off, will provide the "vacuum" Gaussian wave 
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packet with width L = y , where m is the mass of the particle and u the frequency of the oscillator. 
Note that the dispersion time r coincides with the inverse of the frequency of the oscillator. If the 
harmonic oscillator is in the n-th excited state, then a state with (n + 1) "humps" is obtained. 

If we "capture" one of these traveling states at time T\ switching on a harmonic oscillator potential 
with an appropriate frequency wi, it would "freeze" in a harmonic oscillator state (without dispersion), 
until the potential is switched off again at a time T2 and the wave packet is released, traveling again 
as a free wave packet that disperses in time. This way, information might be stored temporally in this 
"oscillator traps" , which can also be used to further manipulating them or even measuring the resulting 
state by means of adequate lasers. The frequency uj\ required to capture the dispersed wave packet might 
be fine tuned in such a way that the wave packet, at the time T\, matches an appropriate eigenstate 
(with the same n) of the harmonic oscillator with frequency uj\ up to a phase (see below). 

Which is the resulting wave function ip(T^) after the whole process, given a state prepared for the 
initial harmonic oscillator ip' '(To)? Obviously, the formal solution will be a product of three evolution 
operators, describing free evolution from To to Ti, harmonic with new frequency u\ from T\ to T2 and 
free from T2 to T3: 



Before proceeding, let us discuss a little bit of notation. In the previous sections, we have used 
lowercase, unprimed letters for quantities referring to the free particle and lowercase, primed ones for 
those of the harmonic oscillator. Now we are indicating physical or true time with capital T's. As the 
Arnold transformation includes a diffeomorphism in time, if we are going to use it as a tool to perform 
calculations, we have to be very careful of not confusing the physical time with that used in the Arnold 
transformations. Recovering the notation of lowercase letters, fix To = to = t' . Then T\ = t% = ^FJ' 
with u\ and U2 satisfying (J8|) at t' Q , and t' 2 = T2, for each evolution to be the initial condition of the next 
one. Let us also call t = T3. 

Regarding the transformations we are going to use, denote by A^^it) the QAT from a harmonic 
oscillator of frequency u performed at time t, and by the unitary time evolution operator for that 
harmonic oscillator. With this notation and the previous choice of classical solutions u\(t') and U2(t'), 
for t = to Aujoito) is the identity. 

Now, the product (|19p can be decomposed by a sequence of QAT's and harmonic oscillator evolution 
operators splitting each free evolution operator as the figure (fl~5l) in Sec. [3] suggests: 



i>{Tz) = U{T Z ,T 2 )U^(T 2 ,T{)U{T U T )1>\T ) . 



(19) 



lf>(t) = ^ 1 ,t a (t)^ 1 (t',t2)^ 1 (t2,tl)^,to(tl)^!,(ti,«d)^'(to) 

= A Wl it2 {t)U' Ul (f, h)!^ (h ) (t[ ,t W (t ) , 



(20) 



where t = , with u\ and U2 satisfying © at t2- This way, in the last expression only evolution 
operators for harmonic oscillators appear. 

The following diagram helps to visualize the setup: 



7 




Here, the solid line indicates the physical process followed and the dashed line the "computational" 
timeline we propose here. The line on the right side denotes the physical timeline. 

An obvious generalization of the proposed method consists in increasing the number of captures and 
releases. On the other hand, depending on the particular initial state, it may be a better choice a path 
of calculation along purely free evolutions and inverse Arnold transformations. 

Note that the evolution will be particularly simple if the frequency uj\ is such that the width of the 
wave packet fits a natural width of this second harmonic potential^. When the frequency of the harmonic 
oscillator is not modified, and the same frequency oj is used to capture the state in the harmonic oscillator 
trap, the resulting state will be a squeezed state (up to a phase) with squeezing parameter r given by 
r = — \og{\5\\) (where 8\ = 1 + iujti), which is negative. This can be seen as a feasible way of producing 
squeezing in trapped states, simply switching off-switching on the trap for a lapse of time t, resulting 
in a squeezing parameter r = — ^log(l + u 2 t 2 ). In fact, a similar way of producing squeezed states in 
Bose-Einstein Condensates was reported in [TS] . 



intuitively, one might say that, in this situation, the particle would be captured in an eigenstate of the oscillator 
Hamiltonian with frequency w\. However, this is not the case: it can be checked that an extra, position-dependent phase 
appears. This fact and the relationship with the "quasistationary" or "pseudostationary" states A. will be analyzed elsewhere. 



8 



5 Summary and comments 

We have reviewed the formulation of the quantum Arnold transformation, particularized it for the 
relationship between the free particle and harmonic oscillator and suggested to use a series of QATs and 
pure harmonic oscillator evolutions to perform the calculations to obtain the resulting wave function of 
a particle that is subjected to the action of turning on an off a harmonic potential. 

It must be emphasized that this method can be easily generalized to different quadratic, even time- 
dependent, potentials and external forces, both changing their time dependence abruptly. Even in this 
case, only classical solutions of each "phase" of the potential and external force are required to construct 
the resulting wave function and evolution operator. We believe that this tool may be useful for analytical 
computations. 
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